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INTRODUCTION 


The effect of suction, through an opening in a surface, upon a stream 
of flowing gas is determined. The quantity of interest is the maximum dis- 
tance, Y, from the surface that the main strean is affected by the applied 
suction. The fraction of this distance Y of the boundary layer thick- 
ness 0 is determined as a function of tubing radius and length, gas pres- 
sure, temperature and stream velocity. The shape of the flow in the high 
velocity streaoa is assumed to be determined by means of a parameter which is 
an unknown function of the distance from the leading edge of the surface. 

The parameter is given values up to and including those that have been shown 
to lead to flow separationC 1 ) . 

MODEL 

The experimental situation is shown in figure 1 where the relevant 
parameters are also defined. The distance Y represents the streanline which 
has a velocity u(Y) in the region prior to the suction. This velocity has a 
zero value at the downstream edge of the suction area. An extensive deriva- 
tion of the relevant equations will not be shown in this report. The de- 
tails are given in reference 1. A few of the pertinent equations will be 
developed as the/ bear on the. proper understanding of the terras in the final 
working equation. 

The starting point is the first order equation for the momentua loss 
thicknessC 1 ) . 

^6^ + ^ ^5 2 + ^ - M^2l" 1,“ “ 0 

ug dx L <^2 J Pj u^2 

where 62 * momentum thickness 
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5^ ■ displacement thickness 
Ug ■ freestream velocity 
T ■ stress at wall 

Hi 

■ mass 

P, ■ density 

0 

Letting u “ 0 and y " 0 at the downstream edge and assuming mass, momen- 
tum and energy are equal to the profile at point 1 and point 2, the follow- 
ing relationships are derived for the quantity of interest, the displacement 
thickness: 


Y 

(5iq) “ I (1- L. ^ 
'' 11 _ 


( 2 ) 

(3) 


where the subscript I stands for the region before suction, II stands for 

the region after suction and Q represents the suction region. Equation 3 

gives a correction to the streamline by a ratio of velocities. This 

correction is called a "suction thickness" as defined in equation 4: 

■Y u 




uS 


dy 


(4) 


The product of the suction thickness and the velocity represents the amount 
of material that can be removed through the opening per unit area. The 
amount of material Q is defined in equation 5 as: 


Q “ v;^ 6q (5) 

In order to perform any of the above integrations a functional form for 
u/u5 is needed. For the purposes of this study, assuming turbulent boun- 
dary layers, the following simplified form is used: 
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iihere o(x) !• the boundary layer thickness. The parameter k is used to 
define the flow characteristi«,» . A value of k above .7 will be beyond 
the point of flow separation. It must be emphasized at this point that one 
of the basic reasons for selecting this functional form is for ease of inte- 
gration and not because of any rigorous justification from fundamental equa- 
tions of flow or empirical evidence. Substituting 6 into 4 and integrating 
\e obtain: 


“ (?) (kTt) 


k + 1 


(7) 


where k ■ k(x) 

Rearranging and solving for Y the following expression is obtained 

1 


Y r ^ ^ 

I - ^ i! 

o L o - 


Eliminating with equation 5 the working equation is obtained: 

1 


( 8 ) 


; ■ [(^) o] 

This equation relates the maximvm distance to parameters that can be readily 
calculated from other equations as will be shown below. Uie only parameter 
that is arbitrarily varied is the shape parameter, k, and that only within 
the limits described above. 

Equation 9 can be generalized by introducing a dimensionless length 
x/L where L is an arbitrary body length. Vhen this is done a new Rey- 
nolds number is defined, 
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pu5L 



(9 a) 


which leads to Che definition of a boundary Chickness 


for the body, 


and Che modified equation 9c. 


Y - 


I k k 

1 k + 1 TTk + 1) 


(k + 1) Q 
u6 



(9b) 


(9c) 


The log(Y) is shown in figure 2 as a function of the velocity ratio 
for various values of with x/L • .5 and k ■ .5. 


EVALUATION OF PARAMETERS 


1. Flow 

Q in equation 9 can be determined using the Poiseuille equation 
for flow through a finite tube modified for end effects. (ref. 2). This 
equation is valid for the conditions used in thio report where the flow is 
viscous and lami.iar. 


Tia^ (?2 - Pi ) 



where 

P2 “ inlet pressure 

Pi ■ outlet pressure 
a * tube radius 

n ■ viscosity 
L » tube length 
p ■ density 


( 10 ) 
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II. Density 


The denalty p is calculated assiming ideal gas behavior. 

P - M.W. (*P/RT) (11) 

where P ■ inlet pressure 

M.W. “ Molecular Weight 
R “ gas constant 
T ■ Temperature (K) 

The assimption of ideal behavior does not present any difficulty since the 
compressibility factor is very nearly equal to one fjf most of the gases of 
concern, in the region of high temperature and low preseure that we are con- 
sidering. 


III. Viscosity 

The temperature dependence of the viscosity is determined from an em- 
pirical equation proposed by Golubev (ref. 3). The viscosity is assumed to 
be Independent of pressure at the low pressures considered here. The Golu- 
bev equation is: 


n 


* 

n T 

c r 

He 

n T 

c r 



965 


Tr < 1 


( .71 + .29/T ) T 
r r 



> 1 


( 12 ) 


1 _ 

T 6 
c 

where T ■ reduced Temperature 
r 

» reduced pressure 
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■ critical temperature 

P ■ critical pressure 
c 


The viscosity at high temperature was also calculated, fot one set of condi- 
tions, using the following relationship proposed by IhodoO) 

-4.058 T 


-.449 T 

- 4.610 - 2.04e ^ 

^ ^ 1/6 ^ - 1/2 „ - 2/3 

C " T M P 

c c 


+ 1.94e 


+ . 1 


( 13 ) 


where the symbols are defined in equation 12. Ihe results of these 
calculations are shown in table 1 for Nitrogen. The two methods are seen to 
be in close agreement with one another and with the experimental values. 

The equation of Golubev was then used in all calculations. 


IV. Boundary Layer Thickness 

The boundary layer thickness is obtained either from shuttle data 
(ref. 4) or from tne following equation: 

o(x) - X / /rIT (14) 

where x is the distance from the leading edge of the surface and Eb* is 
the Reynolds number defined by: 

(15) 

n 

There aro other definitions of the boundary layer for other typ>68 of flow. 
These will not be considered here. The expression used here is the one most 
typically used for flow across a flat plate. 


CALCULATIONS 

In all of the calculations described in this report the flowing gas was 
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afsumed Co b« NITROOEN. Ihe paraneters chac can be adjuaCed are teen to be 
of C«D Cypea. In Che first group are chose variables which affect the 
suction velocity and the second group consists of those that effect the high 
velocity stream characteristics. Ihe suction paraaeters shall be considered 
first. 


A. Suction Paraneters 
1. Variation of Temperature and Pressure 

The effect of temperature on flow velocity is shown in figure 3. The 
decreasing flow velocity, with increasing temperature is due, primarily, to 
the increase in viscosity of the gas. At the highest temperatures the vari- 
ation of the velocity is not as great, attaining a velocity which is rela- 
tively insensitive to further temperature changes. 

In the tunnel experiment, the temperature of the gas undergoes a rapid 
rise of 175*C during the experiment. The observed flow approaches a maximum 
followed by a decrease in flow as the experiment continues. The decrease 
is partially explained, as discussed above, by the increase in the gas 
/iscosity during the experiment. Taking the temperature of the gas to be at 
the midpoint of the temperature rise, the calculated time of arrival at the 
mass spectrometer is 6 seconds, approximately that which is observed. 

The effect of pre*)sure and temperature upon the maximum suction length 
is shown in figure 4. At a given condition of tube radius, length and 
stream velocity increasing pressure and lowering temperature raises the 
likelihood that the material removed from the streon will come from outside 
the boundary layer. The distance Y can be controlled by selecting the 
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appropriate tubing radius and lengthy as shown below. 

2. Effect of Tubing Length and Radius 

In figure 5 the effect of tubing radius and length is sunmarized at one 
particular condition of streamflow. The larger the radius, or the shorter 
the tubing, the higher the suction velocity with a resultant increase in the 
Y/o ratio. An effective sunmary of the tubing parmeters is given in 
figure 6 where Y/o is plotted against a^/L, a measure of the ratio of a?- 
cross-sectional area to tubing length. 

B. Effect of Streamflow Characteristics 

1. Stream velocity 

At high mach flow, the effect of the streamflow velocity upon Y 
diminishes. For a mach number greater than 2 the ratio Y/o does not 
vary appreciably with stremn velocity ae sho%m in figure 7. 

2. Flow Shape 

The flow shape is governed in this model by the parameter k the vari- 
ation of this parameter and its effect upon the Y/o is shown in figure 7. 


SUMMARY 

Wind Tunnel Experiment 

For the wind tunnel data, the effect of the cooling Nitrogen upon the 
boundary layer thickness was ignored. Evaluating equation 8 for the tube 
length (1800 cm), tube radius (.08 cm) and a Mach of 5 we find the Y/o ■ 
.03, well within the calculated boundary layer. Other characteristics of 
the flow, delay time and temperature behavior also agree with calculated 
behavior as discussed in the previous section. This agreement is, quite 



likely, a conaequcncc of the application of the pouaeille equation for tube 
flov than to any intrini ic validity of th«^ aaavaoptiona made abo> t the stream 
flow. 

%utt^e Flight Path 

The shuttle data from reference 3 iras treated with this model. Using 
the boundary layer calculated in reference 3 it was found that with a ori- 
fice radius of .16 and a tubing length of 10 cm. roDOval of laaterial from 
the stream would lie well within the boundary layer calculated in reference 
3. As this boundary layer is quite large it would take an extremely high 
suction velocity to penetrate it. These calculations are not sunmarized 
within this report. 

Evaluation of the Appvoxi' uitions in the htodel 

The greatest amourtt of uncertainty with this approach occurs with the 
assumption foi the velocity is'cio made in equation 6. The basic reason for 
the assumption is the simplification of the calculation. Other velocity 
profiles are being examined. 

The equation, used to calculate density and viscosity are approxima- 
tions. However, at the conditions nf ^.^mperature and pressure considered 
here, the error introduced by them is thought to be small. 

The boundary layer thicki.’ess i*i a difficult quantity to define or cal- 
culate. Other approaches to the calculation do exist and will be tested. 

There seems to be no difficulty with the modified poiseuille equation 
for the rlow in the tube, as long as the flow is viscous. Under conditions 
of '^xcremely low outer pressure the possibility of moleculsr flov occurs. 

The equations will Chen have to bo modified. 
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Table 1. 

Comparisons of calculated 

viscosities. 



Thodos 

Golubev 

liR 

T(K) 

in 

Poise 


1000 

407 

402 

399 

900 

381 

375 

375 

800 

353 

346 

350 

700 

324 

316 

321 

600 

292 

285 

289 

500 

258 

252 

255 

400 

221 

216 

218 

300 

178 

176 

178 
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Figure 1. Assumed flow pattern across a flat plate. 
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Suction velocity (cm/sec) 



T (K) 


Figure 3. Suction velocity as a function of flow temperature. 

Radius => .008 cm, length = 10 cm, Mach // = 5, k = .5, 
Pressure = .2. 
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Length (cm) 


Figure 5. Dependence of suction ratio on tubing length (L) and 

radius (R) . Pressure “ .2, Temperature « 750 K, Distance 
from edge ■ 10 cm, Mach #5, k ■ .5. 
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Mach // 


Figure 7. Variation of suction ratio with stream velocity 

and shape factor. Tube radius ■ .008 cm, Tube length 
■ 10 cm, Temperature « 1000 K, Inlet Pressure ■ .2 atm. 
Distance from edge ■ 10 cm. Outlet Pressure ■ 0 atm. 




